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Abstract 

An acoustic analogy methodology for improving noise predictions in hot round jets is presented. Past approaches 
have often neglected the impact of temperature fluctuations on the predicted sound spectral density, which could be 
significant for heated jets, and this has yielded noticeable acoustic under-predictions in such cases. The governing 
acoustic equations adopted here are a set of linearized, inhomogeneous Euler equations. These equations are 
combined into a single third-order linear wave operator when the base flow is considered as a locally parallel mean 
flow. The remaining second-order fluctuations are regarded as the equivalent sources of sound and are modeled. It is 
shown that the hot jet effect may be introduced primarily through a fluctuating velocity/enthalpy term. Modeling this 
additional source requires specialized inputs from a RANS-based flowfield simulation. The information is supplied 
using an extension to a baseline two-equation turbulence model that predicts total enthalpy variance in addition to 
the standard parameters. Preliminary application of this model to a series of unheated and heated subsonic jets 
shows significant improvement in the acoustic predictions at the 90° observer angle. 

I. Introduction 

This paper discusses application of an improved physics-based methodology in predicting jet noise for heated 
subsonic jets, with model extensibility for future predictions of supersonic heated jets. The goal is to develop a 
methodology that accurately computes the far- field sound for jets at all observer angles and all operating conditions. 
In order to achieve this goal, efforts are ongoing at NASA Glenn Research Center to examine the current theories in 
jet noise and to combine that with the best practices in numerical modeling, all of which must be validated by 
experiments. Physics-based jet noise predictions codes such as JeNo (ref. 1) that depend on their source modeling, 
as well as mean flow and turbulence information from RANS codes, are inevitably going to be sensitive to the 
details of the modeling and the accuracy of their RANS input. JeNo, which is still in its development phase at 
NASA Glenn Research Center, has been validated for unheated axisymmetric jets at subsonic and low supersonic 
conditions. Sample calculations at heated conditions, however, point to a noticeable under-prediction relative to 
measurements. 

To motivate the discussion, consider the 90° jet noise predictions at the Tanna test conditions SP07and SP49 
reference 2, i.e., stagnation temperature ratio of 1.0 and 3.14 (see table II). The predictions (ref. 1) use the so-called 
self-noise source term in the acoustic modeling (eq. (B5), appendix B), and RANS input from the NASA Wind 
code. As shown in figure 1, hot jet spectral predictions tend to be significantly under-predicted with respect to 
acoustic data. Moreover, this observation applies to predictions using CFD inputs that otherwise provide good 
agreement with experimental PIV data for both mean and fluctuating velocity in the developing jet region, where 
most of the noise is produced. For example, a detailed CFD study was conducted utilizing available PIV data 
measured in the developing jet region for a range of jet exit velocities and temperatures. Findings (ref. 3) indicate 
that the turbulent kinetic energy levels and peak centerline locations predicted by a standard high Reynolds number 
variant for the k-e model trend quantitatively well with increased jet temperature. In addition, mean flow mixing 
appears well captured without the need for additional modifications to enhance shear layer mixing growth due to 
density gradients. However, as seen in figure 1, the predicted changes in the turbulence values are insufficient to 
account for the peak noise differences between heated and unheated jets. Consequently, it appears plausible to 
attribute the noise prediction deficit to heat-related sources. 
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Figure 1. — Spectra/data comparisons (dB) using baseline JeNo code for subsonic 
jet cases at 90° inlet angle: (a) SP07 (unheated) and (b) SP49 (heated)(ref. 1 ). 

Noise from heated jets is of particular interest in aeroacoustics due to its far-reaching implications in defining the 
source. Lighthill’s definition of equivalent aerodynamic noise sources Ty = pv,v / + \ p- cjjp j 5,y - x ( y has been a 
subject of wide-spread scrutiny and interpretations. This has led to scaling laws that relate the far-field sound to flow 
parameters such as velocity and temperature in a variety of forms. In the absence of viscous stresses Xy, the 
momentum flux term d 2 pVjVj ld 2 x i Xj is usually considered as the primary source of sound in isothermal jets. Heat 
addition has a multitude of effects on noise that depends on the jet velocity, frequency, and angle. Experimental 
observations (refs. 2, 4, and 5) concur that at high speeds, heat addition results in noise reduction at all frequencies, 
while at low speeds, say U /c x = 0.50 or below, it amplifies the low to mid frequency jet noise with minimal 
effect at the high frequency. In between, say at U jfc m -0.70, the spectrum shows a cross-over relative to the 
unheated case that indicates an enhancement of the low-frequency noise and slight weakening of the high-frequency 
amplitude. 

These observations have compelled a number of researchers to examine the so-called pressure/density term 
pj as a potential additional source and formulate scaling laws that resemble the usual t/law suggested based 
on the first source term. A recent noise data analysis (ref. 6), examines a more general power law AU " where 
amplitude A and power n are parameters that depend on angle and jet temperature. Morfey et al. (ref. 7) revised 
Lush and Fisher’s (ref. 8) version of entropy source term and offered a new two-scaling hypothesis, in favor of if 
and Empower laws, and with some dependence on the mean temperature gradient. Additionally, they prompted to 
account for the mean flow effects using some form of Geometric Acoustic approximation. Lilley (ref. 9) proposed a 
second order wave equation for the pressure fluctuations, similar to Lighthill’s, but replaced with its equivalent 
from the energy equation. His new contribution to the total acoustic power consisted of a Dipole term with if 
dependency, which also multiplied the variance of the total enthalpy. He provided some estimate of the relevance of 
the new source term without actually calculating the required enthalpy fluctuations. 

The predictions provided in the current study are based on a more systematic form of the governing equations 
that separates the background flow from the fluctuations starting from the set of Euler equations. Following an 
analysis similar to Goldstein (ref. 10), and with suitable assumptions about the base flow, these equations are 
combined to form a third-order linear wave operator with pressure fluctuations p ' as its dependent variable, and a 
set of non-linear terms that are second-order in fluctuating variables and are identified as the equivalent sources of 
sound. As seen in the following discussions, the analysis seems to point, among other factors, to a combination of 
fluctuating velocity/enthalpy term as an additional source. 

To improve the RANS-based modeling capability for noise, a set of transport equations is utilized that predict 
the total temperature (or enthalpy) fluctuations in addition to the usual turbulence-related parameters. The approach, 
as discussed in (ref. 11), makes use of a baseline k-e turbulence model demonstrated in past work to predict correct 
mean flow velocity mixing, turbulent kinetic energy, and Reynolds stress variations for heated and unheated 
subsonic round jets. A generalized two-equation scalar variance model is available within this turbulence model 
framework, which utilizes an independent dissipation rate equation as well as the locally available turbulent velocity 
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and time-scale information. Determination of the modeled scalar variance variable is principally achieved by 
selecting the appropriate mean flow quantity gradient for its production source term. In the past, this model has been 
used to predict the local variations in the turbulent Prandtl number and to study its impact on jet thermal mixing 
using an energy variance approach that tracked the variance in static temperature. The current study traces the 
variance of the total enthalpy instead as the parameter of interest for subsequent jet noise predictions. The difference 
between these two parameters becomes increasingly obvious with jet speed, producing a relatively noticeable static 
temperature variance within the jet shear layer even under unheated conditions compared to practically insignificant 
total temperature variance. It is interesting to note that, at least for the jet cases studied, the use of total enthalpy 
variance does not significantly alter predictions for the variable turbulent Prandtl number. The local turbulent 
thermal diffusivity is partially determined from the thermal time-scale; the predicted values for the variance 
dissipation rate are also affected proportionately by the scalar variance production term selection. 

The remainder of this paper is divided into the following sections. The governing acoustic equations and the 
sources of jet noise in the specialized case of a locally parallel mean flow are presented. This builds grounds for 
modifications subsequently introduced into the transport equations that accommodate the flowfield information 
required for noise prediction. A summary of the baseline scalar variance model framework is presented first. Details 
regarding the model development can be found in (refs. 11 and 12). The variances in total and static temperatures 
are then compared for a series of heated and unheated subsonic jets. Comparative results for the subsonic jet cases 
feature significant improvements in predicted jet noise for heated jets. Discussions are currently restricted to the 90° 
observer angle; results for other observer angles will be forthcoming in a future paper. 

Extension of the modeling methodology to supersonic jets is still exploratory. Most flowfield modeling work 
related to jet noise in this area has focused principally on obtaining correct mean flow velocity and/or temperature 
decay rate predictions relative to available measurements. The experimentally observed longer potential core lengths 
for supersonic unheated jets, which decrease as the exit temperature is increased, have prompted use of 
compressibility corrections to the standard turbulent kinetic energy equation for best agreement with data. However, 
a physics-based jet noise prediction methodology additionally requires accurate portrayal of the underlying 
turbulence quantities in the developing jet region, and this remains an area where further research is needed. In the 
present study, we restrict our study of supersonic jets to highlight desired numerical sensitivity results for total 
enthalpy fluctuations in the mixing region. Further efforts for supersonic jets will require studying the impact of 
turbulent compressibility corrections on the predicted jet noise outcome and may likely require use of both 
experimental and numerical (using LES) findings as guidance. 

II. Acoustic Equations 

This section briefly describes the governing acoustic equation, which can best be described as a form of Lilley 
acoustic analogy in a locally parallel mean flow. The details of the source modeling and source/Green’s functions 
convolution integral are given in appendices A and B. 

Assuming that the ideal gas law p = p'-RT holds, each flow variable is divided into its mean and fluctuating 
components 

P=P+P\ P = P + P'> v,. = v,. + v'i , h = h + h ', ^ 

where the over-bar is used to denotes a time-averaged value, and tilde is a Favre-averaged (i.e., mass-averaged) 
quantity, and the stagnation enthalpy is 

h t =h + ^v 2 , h t =h+^v 2 , h' t =h'+^v' 2 . (2) 

The set of Euler’s equations are now linearized using the above definitions, and expressed as two sets of equations 
governing the base flow and its fluctuations or residuals (ref. 10). Additional simplification is achieved only when 
isentropic form of the energy equation is considered and the base flow is described as a locally parallel mean flow 

Vj = 5 /1 f/(x 2 ,x 3 ), p = p(x 2 ,x 3 ), J? = constant (3) 

The governing equations are thus combined into a single third-order wave equation 

Ln’ = T, n' = (4) 

y p 
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where the linear operator L is 
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In view of the ideal gas law, we require 


_ ~ y — — — 2 

P h = -p, y p = pc 

Y-l 


(8) 


Three pairs of brackets in eq. (7) designate three distinct sources. The former terms in each bracket are identified as 
self, shear, and enthalpy sources respectively. Self and shear source terms are second-order in fluctuating velocity, 
and are each complemented by a second term involving the mean flow enthalpy which is not zero even when the 
flow is isothermal, but may be considered small relative to the former. For example, sample predictions of (ref. 13) 


conclude that ^—D 2 

dx .- 


v lVj 


is small relative to other terms in unheated jets. Source terms that consist of the mean 


density gradient are negligible in unheated jets, but could potentially become important as jets get hot. A detailed 
comparison of the individual source terms of eq. (7) will be presented in an upcoming paper for all set points shown 
in table II. It turns out that the major difference between the unheated and heated jets is due to the fluctuations in 
stagnation enthalpy that appears in the third bracket above. This term is first order in velocity fluctuations and its 
radiated sound power dominates that of density gradient term as well as the self-noise source component as jets 
become hot. 

Spectral prediction presented in section VII display the two relatively more significant source terms denoted as 
self- and enthalpy-noise, and are based on eqs. (B5) and (B6) discussed in the appendix (see appendix A for source 
modeling and appendix B for the derivation of the spectral functions). 


III. Enthalpy Variance and Dissipation Rate Model 

The high-Reynolds-number form of the equations governing the transport of the total enthalpy variance and its 
dissipation rate is based in concept on the work of Nagano and Kim (ref. 12) for generalized turbulent thermal 
diffusivity 
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Here, k and s denote the turbulent kinetic energy (TKE) and dissipation rate obtained from a companion two- 
equation turbulence model, a and a, are the molecular and eddy diffusivities for heat, and P k is the production term 
used in the TKE equation. Values utilized for the modeling constants in the variance dissipation rate equation are 
based on results described in reference 1 1 and are listed in table I. 


TABLE I.— ENERGY VARIANCE DISSIPATION RATE COEFFICIENT CONSTANTS 


Cdi 

C d2 

c d3 

c d4 

Gkt 

Get 

2.0 

0.72 

2.2 

0.8 

1.5 

1.5 


Physical interpretation of the quantity k , tracked within the CFD simulation is based on selection of the flow 
variable gradient used for the turbulent variance production term P k For example, Nagano and Kim originally 

solved directly for the static temperature variance. Recent numerical studies modeling helium jets and afterburning 
rockets plumes indicate improved model generalization using an internal energy variance approach. As shown 
below, further extension of the framework to model the total enthalpy variance was found necessary to improve jet 
noise prediction. 

Two modeling approaches were studied for use in hot jet noise prediction. The first utilized the internal energy 
variance as described in reference 11. The second modeled the total enthalpy variance, with appropriate variable 
substitution for the variance dissipation rate equation. The production term is then rewritten for the selected method 
as follows: 


_ de de 
= 2a, , 

dx j dx i 

(internal energy variance) 

(11a) 

dh, dh. 

= 2a, — — , 

dx ; dx: 

(total enthalpy variance) 

(lib) 


Where e is the specific internal energy and h, is the specific total enthalpy. For a calorically and thermally perfect 
gas, the specific heat constants are invariant with temperature, and so the temperature variance can be expressed as 
follows: 


= e'e' = c v 2 (rr), 

(internal energy variance) 

(12a) 

= W,=c p 2 (7777) , 

(total enthalpy variance) 

(12b) 


and c p = y c v . Values for the coefficient constants in the variance dissipation rate equation for both approaches were 
unchanged. 

When extending the variance model to predict thermal eddy diffusivity, an expression for the thermal turbulent 
time scale is needed. Based on results from past work, the following expression was used: 


{ k 7 

'V 

, 

(C } 

'“'Ll 

fk) 

— max 


1 

UJ 

~ 

UJ 

_v e tj 


U J 


C p = 0.09 


(13) 


The thermal eddy diffusivity, applied to the turbulent diffusion term in the mean energy equation and the variance 
production term, is written for model closure as 

a, = C x pkT, . (14) 

As recommended by Chidambaram et al. (ref. 14), value of C r = 0. 14 was used. The local turbulent Prandtl number 
is then defined as the ratio of the eddy viscosity to the thermal eddy diffusivity 

P F = IV /a, . (15) 
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IV. Subsonic Round Jet Flowfield Simulations 

Both the energy variance and total enthalpy variance models were applied to the simulations of several subsonic 
round jet cases, with nozzle exhaust conditions defined by set points within the Tanna matrix. A summary of these 
cases is listed in table II. In past work, these cases were studied extensively to better understand and predict 
turbulence quantities in the developing jet region, where most of the jet noise is produced. Comparisons were made 
with available PIV data collected at the SHJAR facility (ref. 4) at NASA Glenn Research Center. 

The jet flowfields were computed using an axisymmetric 316x171 computational grid, which included the 
internal nozzle region well upstream of the exit plane. The computational domain extended 100 jet exit radii in the 
streamwise direction and 25 jet exit radii radially, sufficiently distant to minimize boundary condition placement 
impact on the shear layer entrainment path-lines. 

Solution convergence was determined for all cases when the L 2 of the residual for the eight transport equations 
dropped five orders of magnitude relative to the initial solution. Solution grid resolution insensitivity was verified 
using a computational mesh with twice the radial grid density and then comparing predicted centerline values 
between the coarse and fine grid solutions. For all cases, the comparisons were identical. 

Nozzle flow boundary conditions were prescribed uniformly as inflow stagnation conditions for pressure and 
temperature imposed well upstream of the nozzle exit. For solution stability, values for turbulent kinetic energy, 
dissipation rate, and scalar variance were monitored throughout the flow domain for numerical exceptions; in such 
cases background ‘‘zero” levels (le-16) were prescribed. 

The boundary layer region is an important aspect of the nozzle flowfield and can impact the strength of the initial 
jet shear layer. However, accurately resolving the turbulence and mean flow strain characteristics through the 
laminar sublayer was beyond scope of the present effort. As a reasonable approximation, a compressible wall 
function formulation was used to analytically surface boundary conditions to points in the inertial region. Grid 
resolution is still required to resolve the momentum deficit of the boundary layer; approximately ten mesh points 
were provided within the boundary layer region near the nozzle exit. Additionally, all wall surfaces were assumed to 
be adiabatic. 

Results for predicting the mean velocity and turbulent kinetic energy using the turbulence modeling framework 
under current study is discussed in detail in reference 3 for the subsonic cases listed in table II. Inclusion of the 
turbulent variance equations had no impact on the prediction of these quantities. 


TABLE II.— TANNA MATRIX SET POINT CONDITIONS 


S.P. 

NPR 

T,JT m 

T/T* 

M a ; 

Mj 

03 

1.197 

1.000 

0.950 

0.50 

0.51 

23 

1.103 

1.810 

1.760 

0.50 

0.37 

42 

1.066 

2.750 

2.700 

0.50 

0.30 

07 

1.861 

1.000 

0.840 

0.90 

0.98 

27 

1.361 

1.922 

1.760 

0.90 

0.68 

46 

1.225 

2.861 

2.700 

0.90 

0.54 

29 

1.900 

2.114 

1.760 

1.33 

1.00 

49 

1.692 

3.138 

2.700 

1.48 

0.90 


*M ac = U j /c. 


V. Jet Scalar Variance Flowfield Results 

Comparison of the resulting static and temperature fluctuations using respectively the energy and total 
temperature variants of the scalar variance equations are shown in figures 2 and 3 As shown, at a given exit velocity 
(i.e., constant acoustic Mach number M ac ). the peak values for both the total and static temperature fluctuations 
increase with exit temperature. For the unheated jets (such as SP03 and SP07), static temperature fluctuations 
increase with jet exit velocity. The variance production term is influenced only by the magnitude of the temperature 
gradient. A drop in exit static temperature with increasing velocity results in larger density gradient, which in turn 
enhances the peak temperature variance. Use of a total temperature variance equation remedies this issue since the 
SP03 and SP07 cases are isoenergetic. It is worth noting that the peak total enthalpy also grows slowly with 
increased jet exit velocity due to turbulent viscous dissipation effects on the total enthalpy (see ref. 15). For 
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example, for a turbulent round jet, the energy equation can be rewritten in terms of the total enthalpy equation with a 
source term proportional to the local velocity gradient as 



For supersonic jet cases, differences in peak magnitudes for the static and total temperature variance are even 
more significant. Table III summarizes the series of unheated supersonic jets simulated during this study. 
Computational grids and nozzle geometries were supplied by NASA Glenn personnel for this effort. As shown in 
the comparisons of figure 4, the peak centerline static temperature fluctuation can exceed 5% of the exit temperature 
for the Mach 1.8 case, whereas the peak total temperature fluctuations remain small, which is desired in jet noise 
predictions within the current framework. The computational assessment studies support the use of total temperature 
variance from RANS flowfield simulations for use in improving jet noise predictions in heated jets. 


TABLE III.— FLOW CONDITIONS FOR 
SUPERSONIC UNHEATED JET STUDIES 
(T a =294A4K) 


Mj 

NPR 

TIT 

j y J oo 

1.185 

2.363 

0.781 

1.400 

3.162 

0.718 

1.500 

3.647 

0.690 

1.660 

4.618 

0.645 

1.800 

5.708 

0.607 




X/Rj 


(a). M ac =0.50 (b). M ac = 0.90 
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(c). M ac =1.33 and M ac =\A% 


Figure 2. — Predicted temperature fluctuations using energy variance model. 
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(c) M ac =l.33 and M ac =\ .48 

Figure 3. — Predicted temperature fluctuations using total enthalpy variance model. 




(a) Energy variance model (b) Total enthalpy variance model 

Figure 4. — Predicted temperature fluctuations for unheated supersonic jets. 
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Figure 5. — Predicted centerline temperature decay 
comparison for SP27 using variable turbulent Prandtl 
number model. 



Figure 6. — Comparison of predicted thermal eddy time- 
scales using total enthalpy vs. energy variance 
model: SP27 case. 


VI. Thermal Mixing Using Total Enthalpy Fluctuations 

The baseline turbulent Prandtl number model utilizes local temperature variance (later extended to energy 
variance) and its dissipation rate to calculate turbulent thermal diffusivity. To determine the impact on using the 
total enthalpy variance for thermal mixing predictions, the SP27 low-speed heat jet was simulated using both 
modeling approaches. For comparison, the jet was also simulated using a constant turbulent Prandtl number of 0.7. 
As shown in figure 5, the centerline thermal mixing as predicted by the total enthalpy fluctuations is similar to that 
predicted by the baseline formulation. For both cases, the predicted total kinetic energy is the same since use of the 
variance model does not impact this quantity. Figure 6 compares the predicted thermal time scales using both 
variance models and shows similar agreement. Based on eq. (13), the predicted turbulent thermal diffusivity should 
therefore be similar. 


VII. Jet Noise Predictions 

Spectral predictions presented here are for 0 = 90° observer angle and are based on equations (B5) and (B6) 
discussed in the appendix. Numerical results are shown for all eight subsonic test points described in table II, and are 
compared with the lossless narrow-band data of reference 4 at a distance of RIDj = 100. In addition to the standard 
flow variables, which include the turbulent kinetic energy and time- and length-scales, this modified version of the 
flow solver permitted input of the local total temperature variance as predicted by the modified scalar variance 
model. Figure 7 shows that at the unheated conditions of SP03 and SP07 self-noise component dominates. With the 
addition of heat, the enthalpy-related noise overwhelms the self-noise at the static temperature ratio of 1.76 and 
2.70; and increasingly more so at the higher temperature. As is pointed out in appendix B, outside the usual 
calibration constants C e and C T that are related to length- and time-scales, the final calibration of eqs. (B5) and 
(B6) requires a pair of constant A and B. The predictions appear quite encouraging and show a consistent level for 
power spectral density at all test conditions examined here. 

As a cautionary remark, due to the differences in the RANS predictions between various flow solvers, 
specifically in quantities such as turbulent kinetic energy and its dissipation rate, and probably the variance in total 
temperature, it in highly unlikely to define a universal set of constants for the above four calibration parameters. The 
final calibration thus becomes dependent on the specific RANS solver. 
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VIII. Conclusions 

An acoustic analogy methodology has been presented in order to improve the far-field noise predictions in 
subsonic heated jets. The governing acoustic equation is an inhomogeneous Pridmore-Brown equation that is 
constructed from a linearized set of Euler equations when the base flow is represented as a locally parallel mean 
flow. The resulting nonlinear terms, all second-order in fluctuating variables, are identified as the equivalent sources 
of sound and are modeled. 

The effect of heat on noise is attributed primarily to a fluctuating velocity/enthalpy source term. In the absence 
of any experimental measurements, a number of simplifying assumptions are required to model this new source. A 
specialized RANS solver supplies the variance in total enthalpy in addition to the standard parameters. This flow 
solver uses a pair of transport equations, governing the variance in total enthalpy and its dissipation rate, coupled 
with the standard two-equation k-t: model. 

Sample 90° spectral computations demonstrate that the boost in source strength is sufficient to correct the earlier 
deficits predicted in the absence of the enthalpy source term. Outside the two calibration parameters related to time- 
and length-scales of turbulence velocity fluctuations, the final calibration is achieved when two additional constants 
are introduced that amount to using a linear combination of the two sources. This is required partially due to the 
neglect of the cross-correlation between various source terms, as well as differences in scales between velocity and 
thermal fluctuations. 

Extension of the above study to supersonic jets is currently underway. Preliminary RANS solutions indicate 
some level of sensitivity in the predicted turbulent kinetic energy with/without compressibility corrections. These 
corrections are normally introduced to predict a better mean flow, and their impact on turbulence requires further 
consideration. 
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Appendix A — Source Model 


The far-field sound spiral density per unit volume of turbulence at y is 


+00 

P 2 (x, y, ©) = | G(x, y, ro)| 2 J J q(y, x)e ,(ax e~ ,k *=dxd\ , 
| —00 


(Al) 


where q denotes a two-point space-time correlation between fluctuating variables at source points A and B separated 
by space £ and time i , and G is the relevant Green’s function (ref. 16). For simplicity, we consider the turbulence 
as isotropic as described by Batchelor (ref. 17), and describe a two-point space-time correlation of the velocity 
components (notation «,■ is now used in place of v/ as a fluctuating velocity component.) 


( u ,)a( u i)b = u 


2 J 1 2 t, 


(A2) 


dZ, 

Consider a fourth-order time-delayed space-time correlation per unit volume of turbulence at y 

R ijke (y,t^) = ^iUju' k K ■ 


(A3) 


The axial component of the above tensor is 


^liii = [u\ j H^x), 

and its spectral density is the four-dimensional wave-number frequency spectrum function 

Am (y,a s ) = J J Am (ulu) e mx e~ ,kA dxdl. 

I r 


(A4) 


(A5) 


The above integral was written in a frame moving with source convection velocity U c , and ofis the source frequency in 
that frame and is related to the observer frequency through the Doppler factor. If turbulence correlation coefficients are 

quasi normal, then the axial component of the above tensor becomes u x u x u x u x = 2 iu x u x | and subsequently 


mu i 


M= 2 ("i) 2 n 

{ f+W , )~^ L r 

V / \ J J J 

\ 1 

A 2 j 2 ^ 


i(M s x -ik.i 


^dxd\ . 


(A6) 


Similarly, consider the two point space time correlation for the enthalpy/velocity source term 


J ij 1 


V " J A 


V 7 h 


t] 


(y,^\\zs ij (yXx)e i ^e- fk ^dxdl. 


I r 


(A7) 

(A8) 


Using Millionshtchikov’s hypothesis (ref. 17) and noting that the autocorrelation of ^ u t h , j is zero in isotropic 
turbulence, 
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( u > h ’t) A (ujh;) B = (u,) A (uj) b ( hf ) A (h;) B . 


(A9) 


Now if the space/time functions governing the enthalpy fluctuations decay in a similar fashion as that of the velocity 
fluctuations then 



-|2 

e-^e^dxdl, 

2 £, 


(A10) 


or 


Following (ref. 16) we assume 



i j_ wf 
V p 


Am • 


/(q,x) = exp 



(All) 


(A 12) 


then 


/ im =^LA 2 H(u s) )N(ke); = J%xp -2f— | +ia> s x dx = ^ -. (A13) 

5ir ' v T o ) l + (a/x 0 /2 ) 


where N{kd.) is a non-compactness factor (ref. 18). 

In general, if the thermal eddy time-scale Xi is different from the turbulence time-scale x„ then the following 
Fourier transform 

H(g> s ) = f exp -j — 1-|— |+/(o‘ s x dx = x* = 2t ° Xi (A14) 

J -°° 1 T * T i J l + («,V/2) T « +T i 


is used in the velocity/enthalpy spectrum function (A1 1). 
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Appendix B — Source Green’s Function Convolution 


A stationary point source with frequency co and location X s (superscript s denotes a source location) is 
considered in defining the Green’s Function (GF) 

L(<ae- mt ) = cle- mt d(x-x s ), 

and 



1 /£(./?— jtf cos0 

e v 

4ttco^ 


CO 

'Y f m f r* , k, 9 j cos m ^cp — cp s 

m = 0 


)• 


(Bl) 


where f"‘(r s ,k, 9) is a solution to second-order compressible Rayleigh operator (ref. 16), corresponding to mode m 
(m is not a power), wave number k = co /c x (to avoid confusion, notation k will be used to denotes TKE in this 
appendix) , and observer angle 9. Superscript s denotes a source location. 

The relevant GF associated with a moving singularity with source frequency of and convection velocity iU c is 



Some level of simplification is achieved when we consider components of turbulence fluctuations in the direction of 
the observer 9 as suggested by Proudman (ref. 19). The above sources of interest would now be designated as 

2 

d u„u n 2 d f h* ") 

D — - and D — u -4- where r| is the spatial separation of the correlation in the direction of the observer 

or|or| dv\\ ' h ) 

(see sketch). 

k co/ Cqq 
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The required spatial derivatives of the Green’s function are 


6^ = cos 06 [ + sin 9 <3 r 
6 ^ = cos 2 0Gj i +sin 2 06 , T + 2sin0cos06 1; . 


where 6 t = -ik cos 0 G and the radial derivative Q applies to the function f m (r,k,Q) . 

The following simplifications are introduced into the radial derivatives using the high frequency solution to the GF 

r; = -ikor , f” = -(kQf r , (B3) 

where 

Q 2 =-2-(l-M ? cos0) 2 -cos 2 0. (B4) 

Poo 

A more accurate representation of function Q should account for mode number 

Qm =|~(l _M ' 5 COS0) 2 -COS 2 0-j^j . 

Since the zero* mode usually makes the major contribution to the GF, the last term on the RHS was neglected in 
approximating the derivatives of the GF. 

Subsequently, the following expression are derived for the far-field spectral density at 90° per unit volume ring- 
source at y 


and 


P 2 se i f ( x ,y,(o) = 


f 1 f 
U nR I 


\Pooy 


P 2 enthalpy (*>T>®) = 


i Y 


4tt R J 


p ' 1 (p 

m = 0 

CO 

(c 2 p 2 E nn )k 4 X(i+5 m0 ).r.rY 


m = 0 


i (^ <D *) = 7 




2k 


2 „ 2 h 


2 nnin’ i 


nnnn 


(t,^) = k 7/2 - 


-(ro"x 0 /2) 


-N(k£). 


(B5) 

(B6) 

(B7) 

(B8) 


The turbulence length- and time-scales are calculated the usual way from the turbulent kinetic energy and its 
dissipation rate as £ = Qk 1,5 /e and t 0 = C T k/e . The final calibration is achieved when the above noise components 
(B5) and (B6) are also multiplied by constant parameters A and B respectively. 
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